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Abstract 

$h . We present an interesting result regarding the implication of truncating the wavepacket 

of the harmonic oscillator. We show that disregarding the non-significant tails of a function 
which is the superposition of eigenfunctions of the harmonic oscillator has a remarkable 
consequence. Namely, there exit infinitely many different superpositions giving rise to the 
same function on the interval. Uniqueness, in the case of a wavepacket, is restored by a 
postulate of quantum mechanics. 
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We analyse the effect of truncating the wavepacket of the harmonic oscillator in the light 
of the frame theory. Such a theory, developed in 1952 by Duffin and Shaffer in the context of 



harmonic analysis [1] has been applied, for over fifteen years, to construct coherent states. More 

£3 . 

specifically, we should mention affine coherent states, also called wavelets, and Weyl-Heisenberg 
coherent states, also known as Gabor frames [2-1 1] . We recall in the next paragraph the general 
definition of frames and a few properties, which is all what we need for the purpose of the present 
effort. For a complete treatment of frames we refer to [12-14]. 

Given a Hilbert space 7i, a family {0 n } n£ N m W is called a frame for 7i if for every / G 7i 
there exists a pair of constants < A < B < oo such that 

A(fJ}<J2\(<l>nJ}\ 2 <B(fJ). (1) 

neN 

The constants A and B are called the frame bounds and the frame condition. From its 
definition it is clear that a frame is a complete set, since the relations (<p n , f) — 0, n G N imply 
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that / = 0. The removal of an element from a frame leaves either a frame or an incomplete 
set. A frame that ceases to be complete if an arbitrary element 0„ is removed, is called exact. 
This last property implies that only exact frames are bases; in the general case the family 
{</>n}neN may be over complete. When the condition A = B holds, the frame is said to be a 
tight frame. Assuming that all the elements n , n G N are normalised to unity, a tight frame is 
an orthogonal basis if and only if A = B = 1. Notice that this implies that a tight frame with 
frame bounds A = B = 1 is an orthogonal basis only if the elements are normalised to unity, 
otherwise it is a redundant frame. 

In this paper we introduce a class of redundant tight frames which are trivially obtained 
by redefining the functions of an orthonormal basis to be zero outside an interval. Such a 
restriction allows us to use the truncated functions to represent a given function vanishing 
outside the identical interval. This has a remarkable consequence, namely, the coefficients of 
the corresponding linear span are not unique. We discuss this important consequence in relation 
to the truncation of the wavepacket of an harmonic oscillator. We show that, by disregarding 
the non-significant tails of a function which is the superposition of the harmonic oscillator 
eigenfunctions one creates a null space. As a consequence, the restriction of the wavepacket to 
the reduced interval can be realized by infinitely many different coefficients giving rise to the 
same function on the interval. This is, certainly, a striking result. Nevertheless, uniqueness can 
be restored by means of one of the postulates of quantum mechanics. 

The paper is organised as follows: In Section II we introduce the construction of tight 
frames for the Hilbert space of functions vanishing outside an interval by simple truncation 
of orthonormal basis functions. In Section III we apply these results to analyse the effect of 
truncating the wavepacket of the harmonic oscillator. The conclusions are drawn in Section IV. 
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2 Building tight frames from orthonormal bases 

The next proposition shows that, for the space of square integrable functions vanishing outside 
an interval, one can construct tight frames by simple restriction of orthonormal functions to 
the corresponding interval. 

Proposition 1. Let {ip n }n&i be an orthonormal basis for L 2 [L\, L2}, i.e., 

(^mAn) = Ip* m (x)l/J n (x) dx = 6 m>n (2) 



and X[i' l 'j the characteristic function for the interval [L' 1; L 2 ], i.e, 

, x) = f 1 ifxeiL'^L^} 
^[L lt L 2 }\ y q otherwise. 

Functions {^' n } n &h obtained as ip' n (x) = X\l' i'j^V'n^), n G with [L^, L 2 ] C [Li,L 2 ], 
constitute a tight frame for the subspace W 0/ square integrable functions vanishing outside 
[L[,L' 2 ]. 

Proof. On the one hand, since {ip n } n eH is an orthonormal basis of L 2 [Li,L 2 ], for all / G 
L 2 [Li, L 2 ] we have 

n/ir=</,/>=E</'^><^/>- (3) 

nGN 

On the other hand, for all / G W it is true that X\l' l']/ = / anc ^ we further have: 



nSN nSN 



which proves that {^ n } ng N is a redundant tight frame for W, since the elements if)' n are not 
normalised to unity. □ 

We prove next that, through the finite subset of frames elements if)[, . . . ,if)' N constructed 
as indicated in Proposition we can construct the orthogonal projection of / onto S = 
span{^i, . . . ,ip N }, restricted to [L' 1; L 2 \. 
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Proposition 2. Let . . . , ip' N be as defined in Proposition^ and for each f G W let us define 
a function f N as: 

N 

f N = Y,^Mf). (4) 

1=1 

The function given in (J3J) satisfies: 



f N = { ^ s ^( x ) $ x e t L l' L 2] 
\ otherwise, 



where Ps stands for the orthogonal projector operator onto S. 

Proof. Since {^i,...,^} is an orthonormal set for L 2 [Li, L2], the orthogonal projection of 
/ G L 2 [Li, L 2 ] onto S is given as 

N 

P s f = YsM^,f)- (5) 



For / G W it holds that 



Then 



i=i 



N N 



i=l i=l 



N 



P s f = J2^MJ), (6) 



i=l 



and the proof follows by multiplication of both sides of the equation by X[l' i!\ ^ 

A convenient property of tight frames is that the coefficients of a linear span in a tight 
frame superposition are obtained by inner products with the frame functions. In our context 
this implies that the coefficients of an orthonormal expansion and the ones to span a function 
in W by means of the frame in Proposition ^ are computed in an equivalent manner. The 
essential difference is that, as discussed below, the coefficients in the frame superposition are 
not unique. 

Let us consider the frame constructed in Proposition^and let us define S' = spanjf/4, . . . , ip' N }. 
For every / G S' C W we have: 

N 



f = J^c n <, with c n =(tfj' n J}. (7) 



n=l 



Since / £ W it is true that c n = (ip' n , f) = (ip n , f) and the equivalence with the orthonormal 
case follows. However, the redundancy of the frame implies that the coefficients in ((7|) are 
not unique. Indeed, since a redundant frame is linearly dependent, the following situation can 
occur: 

N N 

o = E c ^ for Ei<J 2 ^ - ( g ) 

n=l n=l 

Taking inner product both sides of the equation on the left with each if)' m we have 

JV 

= 52<WnM> (9) 

n=l 

that we can recast as 

= Gd (10) 

where G is a matrix of elements g m>n = (ip' m , ip' n ) ,n,m = 1, . . . , N and d a vector, the com- 
ponent of which are the coefficients d n , n = 1, . . . , N. Equations (jHJ) and (fTU|) imply that the 
general form for (JZj) is 

N N 

f = J2^'n + J2 C 'n^ (11) 

n=l n=l 

with c n , n = 1, . . . , N given in (JJJ) and d n , n = 1, . . . , N the components of a vector d G null(G). 

It is appropriate to stress the significance of Proposition ^ when the interval [Li,/^] is 
actually the whole real line. Then, for numerical calculations one is obliged to work on a finite 
domain. An important consequence of this fact will be discussed in the next section. 

3 On the truncation of the harmonic Oscillator wavepacket 

We show here that the results of the previous section are relevant to the analysis of the trun- 
cation of the wavepacket of the harmonic oscillator. To this end we simulate two different 
situations, which are specially devised to illustrate the phenomenon we wish to discuss. 

Let us consider that a normalised to unity function G Ti. is generated as linear super- 
position of the harmonic oscillator eigenfunctions, i.e., 



5"^ 

where we have written ij) n (x) = % Hn = in terms of the Hermite polynomial H n (x) . The 

^/2 n - 1 (n-l)!0r 

coefficients c„, n = 1, . . . , N in (|T2*|) are simulated according to the equation: 

p -0.0032(n-80) 2 

(13) 



'Y^ N i e~°- 0064 ( n ~ 80 ) 2 

We consider N = 160. The left graph of Figure depicts the corresponding function for 
x G [—1,30]. Since the values \fr(x) are very small outside this interval (\l/(— 1) = 8.7 x 10~ n 
and \l/(30) = 2.7 x 10~ 92 ) one does not commit a significant error by calculating expectation 
values using this domain. In order to illustrate this fact let us calculate x and x 2 as 

/40 
\m(x)\ 2 xdx = 12.56967570231863 
■40 

/30 
\^(x)\ 2 xdx = 12.56967570231863 (14) 

/40 
\V(x)\ 2 x 2 dx = 158.4912423027778 
40 

/30 
\m(x)\ 2 x 2 dx = 158.4912423027777. (15) 

The difference in the values of x cannot be observed in the given format and the values of 
x 2 differ only in the last of the 16 digits. Thus one could conclude that, for the purpose 
of computing expectation values, neglecting the tails of the distribution outside the interval 
x G [—1,30] is not harmful. However, the assumption that = for x ^ [—1,30] has a 

tremendous consequence: the coefficients of the superposition (fT2|) are thereby not unique. In 
order to illustrate this we compute the vectors in the null space of matrix G of elements 



.--so 



Qmn = I H m (x)H n (x) — - dx. (16) 

y ' 7-i 1 ; 1 V 2m_1 (^-l)!2 n - 1 («-l)^ 

We use just one of the eigenvectors spanning null(G), say the vector c', to construct the coeffi- 
cients d' n — c n + c? n , n — 1, . . . , 160, with c n as in (|T3*j) and c' n , n = 1, . . . , 160 the components of 
d G null(G). The right graph in Figure [T] plots the coefficients d' n . We use now these coefficients 
to construct the function ty"(x). By calculating 



r-30 

\^"-^\\ = ^lj \tyt( x ) - q( x )\ 2 dx = 1.4198 x 10~ 18 (17) 



we do not see any significant difference in the functions. However \\c" — c|| = ||c'|| = 1, which 
clearly shows that the function are considering can be generated in many different ways on the 
interval [—1,30]. Even considering the interval [—10,40] null(G) is still not empty. It would be 
empty for [—15, 40], though, but only if the maximum number of states N = 160 is maintained 
fixed. 

It should be noted that, although the 'critical' interval depends on the number N and a 
smaller value of N would decrease the length of the critical interval, there is room for ad- 
justments. Indeed, by considering zero the coefficients which do not intervene in the original 
superposition, but allowing the transformation to have larger dimension, the critical interval is 
enlarged. Notice that this opens the possibility of producing the identical function by means 
of states that were not present in the original superposition. The next example illustrates this 
situation. 

Consider that the coefficients of the superposition (|T2|) are now simulated as 

c n = I =, n = l,...,20 (18) 

The corresponding function ty(x) is plotted in the left graph of Figure 121 In this case the 
critical interval yielding lack of uniqueness is included in the main support of ^(x). Hence, the 
restriction to such an interval is not possible. Now, increasing the value of N to 130 for instance, 
and considering c n = 0, n = 31, ... , 130, the function \l/(x) does not change. Nevertheless, 
we have a 130 x 130 matrix G constructed by extending the interval to one containing the 
most significant support of \I/(x) (in this case [—7, 10]). Taking one of the vectors in null(G) 
we construct the coefficients c" by an equivalent process as in the previous example. With 
this coefficients, plotted in the right graph of Figure El we construct the graph on the left of 
Figure El Notice that, the fact that coefficients c" have significant values for n = 31, . . . , 130, 
implies that \l/(x) can be realized by states which were not present in the original superposition. 
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Nevertheless, the average energy is the same. Indeed, 

N N 

e = «, h^) = E E« + o<c ^ + (19) 

?i=l m=l 

and, since by hypothesis X)m=i 4d = 0> ^ follows that 

N N 

e = E E # c ^> = ( 2 °) 

?i=l m=l 

This result stresses the point that was made initially. Namely, that without normalising the 
coefficients, there could be infinitely many different ways of realizing the harmonic oscillator 
wavepacket on a finite interval. However, according to quantum mechanics each |c"| 2 represents 
the probability of finding the harmonic oscillator in the state n. Hence, to be able to maintain 
this interpretation we must impose the normalisation condition on the coefficients c". As a 
consequence, we cannot use a vector d of arbitrary norm. In fact, in order to construct a 
normalised vector we need to consider c" = c + Dc', with d E null(G) and D a constant to be 
determined by the condition ||c"|| 2 = 1. Writing this condition explicitly we have 

||c 7 '|| 2 = ||cl| 2 + D(c,d) + D*(d,c) + \D\ 2 \\d\\ 2 = 1, (21) 

and, since d and c are orthogonal to each other we further have 

|| c "|| 2 = p| 2 + |D| 2 ||c 7 || 2 = 1. (22) 

The value of ||c|| 2 is fixed by the function ^(rc), in our case ||c|| 2 = 1. Thus, the only possible 
solution to the above equation is D = 0. This leads to the conclusion that amongst all the 
possible superposition giving rise to the same function ^{x) on the interval, there is only one 
which is consistent with the physical significance of quantum mechanics. 

4 Conclusions 

An interesting result, arising by limiting the domain of a normalised function which is the 
superposition of the harmonic oscillator eigenfunctions, has been discussed. It was shown that, 
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when restricted to a finite interval, such a function can be realized in many different ways. 
Although the mean values of the physical quantities are not affected in any significant manner, 
they can be the result of infinitely many different combinations of eigenfunctions. Uniqueness 
is restored by endowing the function with the significance of a wavepacket. It was proved that, 
in that case, there is only one set of coefficients that can fulfil the normalisation condition. 




5 10 15 20 25 30 20 40 SO SO 100 120 140 160 



Figure 1: The graph on the left represents the function ^ constructed by using the coefficients given 
in (|13|) and also the coefficients c" plotted in the right graph. 




Figure 2: The graph on the left represents the function ^ constructed by using the coefficients given 
in ()18|) and also the coefficients plotted in the right graph. 
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